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Orbits for satellites of the Global Positioning System ( GPS) have been determined with 
submeter accuracy . Tests used to assess orbit accuracy include orbit comparisons from 
independent data sets , orbit prediction , ground baseline determination , and formal errors. 
One satellite tracked for 8 hours each day shows rms errors below 1 m even when predicted 
more than 3 days outside of a 1-week data arc. Differential tracking of the GPS satellites 
in high Earth orbit provides a powerful relative positioning capability , even when a rela- 
tively small continental U. S. fiducial tracking network is used with less than one-third of 
the full GPS constellation. To demonstrate this capability , baselines of up to 2000 km in 
North America were also determined with the GPS orbits. The 2000-km baselines show 
rms daily repeatability of 0.3 to 2 parts in 10 8 and agree with very -long-baseline inter- 
ferometry (VLBI) solutions at the level of 1.5 parts in 10 8 . This GPS demonstration pro- 
vides an opportunity to test different techniques for high-accuracy orbit determination 
for high Earth orbiters. The best GPS orbit strategies included data arcs of at least 1 week , 
process noise models for tropospheric fluctuations , estimation of GPS solar pressure co- 
efficients , and combined processing of GPS carrier phase and pseudorange data. For data 
arcs of 2 weeks , constrained process noise models for GPS dynamic parameters signifi- 
cantly improved the solutions. 


I. Introduction 

The Global Positioning System (GPS), expected to be fully 
operational by the early 1990s, will consist of 24 satellites 
evenly spaced in six orbit planes at an altitude of about 
20,000 km. Knowledge of GPS orbits will provide the basis 
for highly accurate ground and satellite user positioning. A 
wide variety of users will benefit from this positioning capa- 
bility. GPS will be used at NASA’s Deep Space Network (DSN) 
stations in conjunction with very-long-baseline interferometry 
(VLBI) radiotelescopes for atmospheric calibrations, precise 


ground station position determination, monitoring of Earth 
orientation changes on time scales of less than one day [1] , 
and possibly time synchronization at the nanosecond level. 
Differential GPS-based accuracies for high Earth orbiters are 
expected at the several-meter level for altitudes of 5,000 to 
40,000 km [2] . Spacecraft maneuvering near and docking 
with the Space Station will carry GPS receivers and will use 
GPS signals for real-time and near-real-time navigation and 
guidance. Low Earth orbiting spacecraft such as TOPEX/ 
Poseidon [3] and the Earth Observing System platforms will 
have orbit determination available in post-real time to an 
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accuracy of better than 10 cm with kinematic smoothing 
techniques using advanced receivers to track GPS satellites 
simultaneously with a worldwide GPS ground tracking net- 
work. The relatively low cost and convenience of GPS ground 
receivers have created many new opportunities to monitor 
cm-level crustal motions in geologically active regions. Very 
dense ground networks may achieve accuracies equaling or 
surpassing those available from other generally more restric- 
tive techniques such as VLBI or satellite laser ranging. 

The GPS applications with the most stringent requirements 
include the DSN applications, subdecimeter low Earth orbit 
determination, cm-level measurements of Earth crustal motion, 
and cm-level monitoring of changes in Earth orientation. To 
reach these goals, GPS orbits will have to be determined to 
better than 50-cm accuracy. The Jet Propulsion Laboratory 
has been developing and testing GPS orbit estimation software 
and techniques for several years with the goal of demonstrat- 
ing the capability for high-accuracy orbit determination. GPS 
data from several field experiments in 1985 and 1986 have 
been used to determine precise GPS orbits. Although only 
seven developmental GPS satellites were operational and the 
ground fiducial tracking network was limited to sites in the 
continental United States, covariance studies indicated that 
with currently available GPS receivers and antennas it should 
be possible to produce orbits for well-tracked GPS satellites 
accurate to 1 m. Achieving this 1-m accuracy capability is a 
major milestone on the road to ultrahigh-precision GPS 
applications. 

In this article we present results demonstrating submeter 
accuracy for the GPS orbits determined from these field tests. 
Ground station coordinates were estimated simultaneously 
along with the orbit parameters. Accuracy of better than 3 cm 
has been achieved over baselines up to 2000 km, proving that 
GPS is already a very powerful technique for precise position- 
ing over continental distances. 

II. Data Acquisition and Processing 

A series of GPS field experiments took place March and 
November 1985, June 1986, and January 1988. These experi- 
ments were organized by JPL and were cooperative ventures 
with several different organizations participating. In the 
March 1985 experiment, data were collected for about 1 week 
at ground sites in the continental United States only. In 
November 1985, the tracking network also included three 
sites in Mexico, and the experiment lasted for about 2 weeks. 
The June 1986 experiment covered a 3-week period and in- 
cluded sites in the Caribbean region as well as a dense network 
of stations in Southern California. The data from the experi- 
ments up through 1986 has been processed, with the January 
1988 data expected to be distributed shortly. This article re- 


ports analysis based on the data collected during the November 
1985 and June 1986 experiments. Results from the March 

1985 experiment have been reported earlier [4] , [5] . 

A. The November 1985 and June 1986 Data Sets 

The November 1985 GPS experiment took place from 
November 12 through November 24. The June 1986 experi- 
ment spanned about 3 weeks; the results presented in this 
article are based on the first part of the experiment from June 2 
to June 10. Figure 1 shows the locations of the ground track- 
ing sites, which represent a subset of the total number that 
participated in the experiments. TI 4100 GPS receivers [6] 
were operated at most of the sites. SERIES-X(7) receivers 1 
built at JPL were used at Mojave and at Owens Valley Radio 
Observatory (OVRO). In November 1985, water vapor radiom- 
eters (WVRs) were available for wet tropospheric delay cali- 
brations at OVRO, at the Mexican sites, and at Mojave for part 
of the experiment. In June 1986, WVRs were used at Haystack, 
Mojave, and some of the Caribbean sites. Dry tropospheric 
delay calibrations were computed from surface measurements 
of barometric pressure. For receivers located at sites where 
WVRs were not available, wet tropospheric delay corrections 
were computed from surface meteorology data. Single-day-arc 
baseline and GPS orbit solutions were generated for examina- 
tion of residuals and to check the quality of the data. Novem- 
ber 12 and 17 were excluded because of data outages and other 
difficulties at some of the fiducial sites. For multiday-arc solu- 
tions, the November 1985 data set, covering 12 days, was 
divided first into two arcs of 7 and 5 days (November 13-19 
and November 20-24). This was necessary due to a maneuver 
of more than 100 km which took place on November 20 during 
which Navstar 4 was moved to a new orbit. Eventually, we 
modeled and solved for the maneuver as described below, and 
a single long data arc covering November 13-24 was constructed. 
During these GPS experiments, periods of common ground 
visibility lasted about 6 to 8 hours. Most satellites were visible 
continuously from a given ground station for only about 3 hours, 
so several times during the tracking period the receivers switched 
to a new combination. Navstar 8 was unusual in that it was 
visible for up to 8 hours from most of the ground sites. 

Because of the short tracking periods (relative to the GPS 
orbital period of 12 hours) from a limited network of ground 
sites, orbits determined from single-day passes in the 1985 and 

1986 field tests were significantly weaker than those deter- 
mined from multiday arcs. The additional strength from the 
multiday arcs derives mainly from the visibility of the satel- 


1 R. B. Crow, F. R. Bletzacker, R. J. Najarian,G. H. Purcell, J. I. Statman, 
and J. B. Thomas, “Series-X Final Engineering Report,” JPL Inter- 
nal Document D-1476, Jet Propulsion Laboratory, Pasadena, Cali- 
fornia, 1984. 
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lites from the ground during multiple orbit revolutions. With 
observations over more than one revolution, the orbital periods 
are more accurately determined and therefore the positions of 
the orbital nodes are more precise. The down-track satellite 
components benefit especially from the multiday arcs, as is 
sometimes manifested by a corresponding improvement in 
the eastern component baseline accuracy [5] . A second advan- 
tage of multiday-arc solutions is the \fn improvement in pre- 
cision, where n is the number of measurements, which can 
apply to any orbital or baseline component. As more satellites 
are launched and the tracking network expands geographi- 
cally, shorter arcs will achieve the same level of orbit accuracy. 

B. Data Collection and Processing 

Both GPS carrier phase and pseudorange data were received 
at all the sites equipped with TI receivers. Carrier phase only 
was used from the Series-X receivers, which are codeless. The 
carrier signals at LI and L2 bands (1.227 and 1.575 GHz) are 
modulated by a pseudorandom noise code called the P code, 
which operates at 10.23 MHz. Continuously tracked GPS 
carrier phase provides a very precise measure of range change , 
while the P code provides a measure of absolute range . The 
pseudorange is considerably noisier than the carrier phase data 
type, and in this experiment the pseudorange was corrupted 
by errors due to multipath signals. The GPS observables at the 
two L-band frequencies are linearly combined to remove the 
portion of ionospheric delay which varies as the inverse square 
of the frequency. For more details about the characteristics 
of the GPS signals, see [7] and [8] . Hereafter, the terms 
“carrier phase” and “pseudorange” refer to these linear com- 
binations of LI, L2 and Pi, P2. The GPS data were processed 
with the GPS Inferred Positioning SYstem (GIPSY) orbit 
determination and baseline estimation software, which was 
completed and tested at JPL shortly after the data were 
collected for the 1985 experiments. 

III. Orbit Determination Approach 

The JPL orbit determination software utilizes a pseudoepoch- 
state U-D factorized Kalman filter 2 [9] . The filter works as a 
batch sequential program with the option to model parameters 
as first-order exponentially correlated process noise, also com- 
monly called colored noise. The GIPSY software uses the 
J2000 reference system with observation partials for param- 
eters computed relative to the satellite epoch states [10] . 

For the November 1985 data set, the nominal GPS ephem- 
erides were obtained by using broadcast orbits as initial values 


2 S. C. Wu, W. I. Bertiger, J. S. Border, S. M. Lichten, R. F. Sunseri, 
B. G. Williams, P. J. Wolff, and J. T. Wu, “OASIS Mathematical 
Description, v. 1.0,” JPL Internal Document D-3139, Jet Propulsion 
Laboratory, Pasadena, California, 1986. 


and iterating to improve those orbits (and remove large unde- 
sired offsets resulting from different coordinate frame conven- 
tions) with a small subset of the data. During the iteration, 
GPS solar radiation pressure coefficients were also determined 
and these coefficients and the new orbits were used as nominal 
values for the more comprehensive precise orbit filter runs. 
In June 1986, the hybrid postfit ephemeris from the Naval 
Surface Weapons Center (NSWC) was used for nominal trajec- 
tories and a similar iteration was performed prior to the final 
filter solutions. 

The high-precision orbit determination strategy used with 
the November 1985 and June 1986 GPS data was based on 
previous experience with the March 1985 GPS experiment 
[5] . A key aspect of the orbit estimation process is the fidu- 
cial concept, where three or four receivers with well-known 
coordinates in a consistent reference frame are held fixed 
while all other parameters, including orbital states and coor- 
dinates of the nonfiducial sites, are estimated simultaneously 
in the filter. Reference 5 and the references therein discuss the 
fiducial concept as well as alternative approaches. The fiducial 
receivers for the experiments discussed in this article were 
collocated at VLBI sites. Haystack, Richmond, and Fort Davis 
were generally treated as fiducials. In some runs OVRO was 
held fixed, and in others Hat Creek was used as a fiducial so 
that one of the normally fixed receivers could be estimated 
using the GPS data in order to test the internal consistency. As 
more GPS satellites are launched and the ground tracking net- 
work is expanded from North America to include stations on 
other continents, we expect that fewer fiducial constraints 
will need to be applied and more station location parameters 
will be determined from the GPS data. 

A. Clock and Bias Parameters 

The results presented here are based on estimation of sta- 
tion and GPS clocks as white process noise. At each measure- 
ment epoch, each active clock is assumed to have a value un- 
correlated to its value at other epochs. Although some station 
clocks were running off hydrogen masers and most of the GPS 
clocks have well-characterized behavior typical of rubidium 
and cesium atomic standards, this extra information was not 
used. The white noise clocks were estimated simultaneously 
with the other parameters in the filter. When the process noise 
model for clocks is white noise, the results are virtually the 
same as would be achieved with double differencing [11] for 
clock elimination. 

The GPS carrier phase, when continuously tracked, pro- 
vides a very precise measure of range change from measure- 
ment epoch to measurement epoch. The absolute phase, how- 
ever, and hence the absolute range from transmitter to receiver, 
is ambiguous by an integral number of wavelengths. For each 
station-satellite tracking pass, a carrier phase bias parameter is 
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estimated, along with the other adjusted parameters, ignoring 
the integer constraint on its value. Over a period of hours, the 
signature of the range change precisely measured with the car- 
rier phase enables’the orbit to be determined. The pseudorange, 
on the other hand, provides a more direct range determina- 
tion but is much noisier and more susceptible to multipath 
errors than the carrier phase data type. The ultimate accuracy 
would be reached if carrier range were available. Carrier range 
is a range determined from carrier phase with the bias ambi- 
guities all resolved; it has the best features of both carrier 
phase (low noise, low multipath) and pseudorange (absolute 
range measure and geometric strength). Successful carrier 
phase ambiguity resolution has been reported over baselines of 
up to 2000 km with bias fixing or bias optimizing techniques 
applied to single-day arcs [12], [13]. The results reported 
here were achieved by estimation of the white noise clock and 
carrier phase bias parameters with a very large a priori uncer- 
tainty, so that their solutions were basically unconstrained. Bias 
fixing was not used to resolve the carrier phase ambiguities. 

B. Tropospheric Delay Fluctuations 

The troposphere was modeled as a spherical shell which 
adds a delay along the GPS signal path: 

o *»..*.(•) 0) 

where p z is the zenith tropospheric path delay and R is an ana- 
lytic mapping function [14] to map zenith delays to line-of- 
sight path delays at elevation 0. There are two components to 
the tropospheric delay-the wet and the dry, denoted here 
with subscripts w and d. The dry component can be deter- 
mined under the assumption of hydrostatic equilibrium using 
the ideal gas law for dry air to better than 1 cm [15] . The wet 
delay component, although considerably smaller than the dry, 
exhibits greater time and spatial variation and is much more 
difficult to determine accurately. 

WVRs were operated at some of the GPS tracking sites 
alongside the GPS receivers. WVR calibrations are believed 
accurate to 2 cm or better for determination of the wet zenith 
path delay [15] . The algorithm described in [16] was used to 
determine these calibrations. At the other sites, surface meteo- 
rology (SM) measurements (temperature, air pressure, and rela- 
tive humidity) were used for the wet zenith delay calibration 
[17] . The SM calibration is much less reliable than the WVR 
calibration because the surface meteorological conditions are 
not always well correlated with the total atmospheric water 
vapor content [15] . 

Wet zenith delay corrections to the calibrations (WVR or 
SM) were estimated with the GPS data for all GPS tracking 
sites. For most sites with WVRs, a constant wet zenith delay 


parameter with an a priori constraint of 3 cm was estimated 
for each 8-hour tracking day. For sites using SM calibrations, 
the wet zenith delay was estimated daily with a 20-cm a priori 
constraint. In addition, stochastic residual delays were esti- 
mated for SM sites in order to remove signatures which could 
result from temporal variations in the troposphere, time- 
varying errors in the SM calibrations, errors in the mapping 
function, or spatial inhomogeneities due to azimuthal asym- 
metry in the water vapor content. In some cases, tightly con- 
strained process noise troposphere residual delay parameters 
were also estimated for the WVR sites. 

C. Troposphere Process Noise Models 

The stochastic model in the GIPSY filter is for a first-order 
exponentially correlated process noise [9] . The measurements 
are processed in discrete time segments, known as batches. In 
each batch, process noise parameters are modeled as piecewise 
constants. At the end of a batch, a process noise time update 
adds noise to the covariance matrix and thus causes the time- 
varying behavior of the stochastic parameters. The process 
noise time update for the ;th batch maps the estimates and 
covariance for the stochastic parameters into batch / + 1 : 

P /+ , = M/P/ + w y (2) 

where p y is a vector of estimates for the stochastic parameters 
and M is a diagonal process noise mapping matrix. The process 
noise w ; - is a random process with zero mean and 

E ( w / w k) = Q V 0) 

where Q is the covariance matrix diagonal and 8 ; - k is the 
Kronecker delta function [9] . The diagonal entries of M are 
given by 

m ij = ex p [-(' /+ i -9 lT n ] (4) 

where tj is the start time for the jth batch and r i j is the time 
constant for the ith stochastic parameter at the/th batch. The 
corresponding diagonal entry in the matrix Q is 



where a iss , the steady-state sigma for the ith stochastic param- 
eter, is the noise level that would be reached if the system 
were left undisturbed for a time much greater than r. The 
process noise model for each parameter is fully specified by 
o ss and r, which can also vary with time, although the sub- 
script j has been left off o ss for simplicity of notation. There 
are two special limiting cases: white process noise, and a ran- 
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dom walk. For white process noise, r = 0, m = 0, and, as can 
be seen in Eq. (2), the a priori covariance for the process noise 
parameters, p, is completely reset at the end of each batch, 
including zeroing of off-diagonal terms and inserting q for the 
variance on the diagonal. For white noise, the process at each 
time step is independent and uncorrelated with the process at 
other time steps. The opposite case is the random walk. Here 
both o ss and r are unbounded, since a steady state is never 
reached and r = °°. For the random walk, however, q is still 
defined by Eq. (2), where M is now equal to the identity 
matrix. The Allan variance [18] , a 2 (Af), which is often used 
to characterize clock and atmospheric fluctuations [19], is 
directly related to the random walk q: 

o 2 a (A/) = q/At 2 (random walk) (6) 

A wide range of process noise models has been tested on 
the March and November 1985 and June 1986 GPS data sets 
[5] , [15] , [20] . The random walk zenith tropospheric delay 
models with in the range 2 to 4 X 10~ 7 km-s” 1 / 2 (for SM 
sites) produced the best daily baseline repeatability, agreement 
with VLBI, and orbit repeatability. When only constant zenith 
delay parameters were estimated, orbit and baseline accuracies 
were worse by about a factor of 2. The value of \/q adopted 
for most of the November 1985 and June 1986 analyses was 
2 X 10 -7 km*s~ 1 /' 2 , corresponding to about 6 cm variation 
over a 24-hr period. There is evidence from VLBI residuals 
[19], [21] that tropospheric delay and delay-rate fluctuations 
can be well modeled as random walks for At greater than a few 
hundred seconds. Since the GPS data were compressed to 
300-sec intervals, the use of random walk tropospheric fluc- 
tuation models for GPS is consistent with the VLBI findings. 
As discussed in [5] , however, the GPS process noise param- 
eters were used to estimate a fluctuating residual correction to 
calibrated data (WVR or SM), and it is not clear that this 
quantity will have the same stochastic characteristics as the 
tropospheric fluctuations themselves. It was assumed that the 
spectral characteristics of both the tropospheric fluctuations 
and the residuals after calibration would be similar. 

D. Solar Radiation Pressure and Other 
Nongravitational Forces 

Multiday arcs covering 1 to 2 weeks were used to achieve 
the highest GPS orbit accuracy. With a global tracking system 
equipped with high-performance GPS receivers and a full 
24-satellite GPS constellation, covariance analysis predicts 
GPS orbit accuracy well below 1 m after just several hours of 
tracking [22]. However, as can be seen from Fig. 1, the 
ground tracking network during 1985 and 1986 was limited, 
with fiducials in North America only. The seven GPS satellites 
that were operating at that time by design tend to converge 


over the southwest United States. This further limited com- 
mon ground visibility to less than 8 hours per day and reduced 
the geometrical strength of the system. In addition, the pseu- 
dorange available from the TI receivers that were used in 
1985 and 1986 was highly contaminated by ground multipath, 
thereby raising the effective measurement noise. Because of 
these factors, multiday arcs were necessary to achieve the 
desired improvement in ephemeris accuracy. 

With single-day (8-hour) orbit solutions, there is very little 
sensitivity to errors in the GPS solar radiation pressure coeffi- 
cients [4] . However, for multiday arcs with multiple orbit 
revolutions, the orbital period is much more accurately deter- 
mined and therefore the results become sensitive to down- 
track errors resulting from integration of accelerations due to 
solar radiation forces acting on the GPS satellites, mismodeled 
solar radiation, or unmodeled forces such as thermal radiation 
from the spacecraft body. The GPS Block I ROCK4 model was 
used to represent accelerations resulting from solar radiation 
pressure. As described in [5] and references therein, ROCK4 
models 13 surfaces on the satellites according to their size, 
curvature, reflectivity, specularity, and absorption character- 
istics. The model as implemented in GIPSY allows for adjust- 
ment of three parameters: G x , G yy and G z . G x and G z are 
scaling factors in the local spacecraft jc and z directions, where 
the z axis is positive along the antenna toward the center of 
the Earth, the y axis is along the solar panel support beam 
normal to the spacecraft-Sun-Earth plane, and the x axis com- 
pletes a right-handed coordinate system. G y represents a con- 
stant acceleration in the jy-axis direction, often referred to as 
the y-bias parameter [23] . 

In principle, the G x and G z parameters should have the 
same value if the spacecraft were perfectly aligned and the 
model were correct. GPS orbits have been determined in the 
past with estimation of only G y and one parameter (desig- 
nated here as G xz ) to represent both G x and G z [24] , [25] . 
For the multiday arcs determined with the March 1985 data 
[5], G x , G y , and G z were estimated independently as con- 
stants for each satellite with the intention of adding an extra 
degree of freedom ( G x ^ G z ) to absorb unmodeled accelera- 
tions and known deficiencies in the ROCK4 model, which are 
thought to amount to as much as 4 m error over a 14-day pre- 
diction interval [23] . This strategy was successful for data 
arcs up to about 1 week long. However, for longer arcs of up 
to 2 weeks, a noticeable and systematic degradation in daily 
baseline repeatability occurred with the three-parameter con- 
stant solar coefficient approach. A new approach was adopted 
in which two constant solar pressure coefficients were esti- 
mated, G y and G xz , along with two tightly constrained pro- 
cess noise parameters, G x and G z . With this approach, daily 
baseline repeatability continued to improve as the arc was 
lengthened. 
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As an alternative to the use of process noise on G x and G ZJ 
small colored noise accelerations were introduced for the long 
November 13-24 arcs. These consisted of three parameters 
representing a constant thrust in the directions of down-track, 
cross-track, and altitude. 

Both stochastic solar pressure parameters and stochastic 
thrust parameters gave similar repeatability results for the 
2-week data arc, Fig. 2. The results were insensitive to the time 
constant, r, as long as the value of q was held about constant 
(see Eqs. 4 and 5). For G x and <7 Z , r was varied from 1 to 28 
days with q « 9.3 X 1G~ 6 , and the repeatability results were 
essentially the same. The results for the thrust parameters were 
more sensitive to the time constant. The best results were not 
obtained until the time constant was made less than 2 days 
and the value oiq^lX 10~ 28 km 2 /sec 4 . Similar results were 
obtained when the time constant for the thrust parameters 
was made as short as 0.5 day. A single stochastic thrust in the 
direction of the satellite long track also yielded good repeat- 
ability. With the present data, which consists of relatively 
short arcs at the same time each day, one cannot determine 
whether the effects are due to solar pressure mismodeling or 
other random or systematic accelerations acting on the GPS 
spacecraft. One possibility is that a solar pressure stochastic 
model is physically the correct approach but the stochastic 
fictitious thrust model works as well, as long as the thrust time 
constant is short enough to absorb the daily variations implicit 
in the solar radiation pressure signature. Note that another 
approach to this problem has been proposed in which ficti- 
tious force parameters with 24-hr resonances are estimated to 
remove solar pressure, gravity, and other dynamic errors 
which tend to repeat daily [28] . It is hoped that the sources 
of these forces can be isolated when more data and global 
tracking are available. This may be possible with the CASA 
UNO data set [29] . 

E. GPS 4 Maneuver 

A maneuver was performed on GPS 4 (SV 8) at approxi- 
mately UT 0320 November 20, 1985. In a few hours the space- 
craft state was changed by over 100 km and it was moved to a 
slightly different orbit. This time was in the middle of the 
experiment but outside the data collection interval, which was 
from about UT 1 100-1900 each day. In order to perform con- 
tinuous orbit determination over the entire experiment inter- 
val (November 13-24) for all the satellites simultaneously, we 
had to model and estimate the maneuver. The NSWC param- 
eterized the maneuver as a constant thrust over 5 minutes. We 
used a four-parameter model with three instantaneous velocity 
changes and a time-of-burn parameter, which we refer to as an 
impulsive motor burn. The impulsive motor burn allows for an 
instantaneous change in position and velocity while leaving the 
acceleration unchanged. Let AV be the vector velocity incre- 
ment, with three components AF^, AF C , and AV L , where 


subscripts H, C, andL refer to the local spacecraft coordinates 
of altitude, cross-track, and down-track. Let AT B denote the 
time of burn, e.g., the interval over which the burn is applied. 
Then the vector change in position is 

Ar=i-AVr £ ={aAr| (8) 

The acceleration, a, defined by AV and A T Bi is the equivalent 
constant acceleration that would be experienced over the 
interval A T B . The solution for the maneuver parameters AV H , 
AF C , AF l , and A T B was obtained by iterating over a data arc 
spanning November 18-22. The maneuver solution converged 
after three iterations and showed meter-level agreement with a 
separate solution for the GPS 4 position using only data col- 
lected after the maneuver (November 20-24). This nominal 
representation for the maneuver was further refined in later 
longer-arc high-precision GPS orbit solutions. It is interesting 
to note that the maneuver solution did not converge initially 
unless we constrained G x = G z for GPS 4. 

IV. Assessment of Orbit Accuracy 

Five criteria were used to assess the accuracy of the GPS or- 
bits determined from multiday arcs in November 1 985 and June 
1986: (1) orbit repeatability; (2) orbit prediction; (3) daily 
baseline repeatability; (4) agreement between GPS and VLBI- 
determined baselines; and (5) formal errors from the orbit 
filter. 

A. Orbit Repeatability 

Orbit repeatability indicates the precision and, to some 
extent, the accuracy of the GPS orbits. Figure 3 illustrates 
how orbit repeatability was computed for the November 1985 
experiment. The purpose of orbit repeatability is to compare 
orbits determined independently without any common mea- 
surements and then compute the rms difference over a time 
interval during which no data were used for either of the two 
solutions. Figure 4 shows the mean of the rms computed for 
all seven GPS satellites over a 6-hr interval on November 17. 
Figure 4 also shows the significant improvement attained when 
pseudorange is processed with the carrier phase and stochastic 
tropospheric delay models are used. From the formal errors, it 
appears that the pseudorange contributed little geometric 
strength to the orbit solutions, since the rms scatter of indi- 
vidual measurements was 100 to 300 cm for 6-minute mea- 
surement intervals, due mostly to ground multipath. However, 
when the pseudorange and carrier phase are processed together 
and a common clock is estimated, the pseudorange provides 
a priori knowledge of the clock and carrier phase bias param- 
eters at the several-nanosecond (100-300 cm) level, signifi- 
cantly improving the orbits. 
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GPS 6 and 8, the two satellites with the most data and best 
ground viewing geometry, had formal errors below 1 m for 
most of the November experiment and had significantly lower 
orbit repeatability rms than the other satellites. Figure 5 shows 
the repeatability computed over a 24-hr interval on Novem- 
ber 17 when no measurements were used for either orbit solu- 
tion. Since the two solutions being compared were determined 
from independent data sets (Fig. 3), it is concluded that sub- 
meter orbit precision has been demonstrated for these two 
satellites. 

B. Orbit Prediction 

Orbit prediction is a stringent test of orbit accuracy, since 
the estimated spacecraft position and velocity are mapped out- 
side of the measurement interval to give a predicted satellite 
state. Orbit errors tend to be magnified in this mapping pro- 
cess. The accuracy of the orbit models used for mapping are 
also tested by the orbit prediction test, in addition to the 
accuracy of the satellite ephemerides. Figure 6 shows the pre- 
diction test we applied to our November 1985 multiday-arc 
orbits. The average rms errors for four of the satellites were 
0.7 m, 0.8 m, and 1.7 m in altitude, cross-track, and down- 
track components. GPS 8, which was tracked longer than any 
of the other satellites, had a prediction rms error well below 
1 m, even when mapped more than 3 days into the second arc. 
For the first 6 hours of the prediction interval, the rms error 
was 50 cm or less for all three position components. These 
results are shown in Figs. 7 and 8. 

C. Daily Baseline Repeatability and Agreement 
Between GPS and VLBI Solutions 

To further assess the GPS orbit accuracy determined with 
multiday arcs, we have examined daily baseline repeatability as 
well as agreement with independent VLBI baseline measure- 
ments over continental distances (1000 to 2000 km). For the 
November 1985 experiment, we examined baselines between 
Hat Creek, CA and Fort Davis, TX (1933 km); Mojave, CA and 
Fort Davis, TX (1314 km); and Richmond, FL and Haystack, 
MA (2046 km). Daily repeatability was computed as the rms 
about the weighted mean of the daily baseline solutions deter- 
mined simultaneously with one multiday-arc orbit solution. 
For these baselines in North America with good common visi- 
bility of the GPS, the rms scatter was 0.3 to 2 parts in 10 s of 
baseline length for all vector components. Agreement with 
VLBI was 0.3 to 1.5 parts in 10 s for baselines with the same 
type of GPS antenna at both ends. Figures 9 and 10 show the 
results for the 2000-km baselines. 

For determination of the Hat Creek-Fort Davis baseline, 
Fort Davis, Richmond, and Haystack were held fixed as fidu- 
cial stations. For the Richmond-Haystack baseline solution, 
a separate filter run was made with Hat Creek, Fort Davis, and 


Haystack fixed as the fiducial reference points. Note that in 
the case of the Richmond-Haystack baseline the Hat Creek 
fiducial does not have data in the second half of the 2-week 
data arc, reducing the number of fiducials to two for the 
second half. The data from the first half has already deter- 
mined the GPS positions sufficiently that the lack of a third 
fiducial does not degrade the quality of the solution. Although 
high-quality orbits are a prerequisite for good precision and 
accuracy over long baseline distances, there are other factors 
that can affect baseline accuracy aside from GPS orbits. For 
example, local site vectors between the GPS and VLBI an- 
tennas sometimes are inaccurate by several cm. One such local 
survey error was recently discovered at OVRO, leading to a 
5-cm discrepancy between the GPS and VLBI Mojave-OVRO 
baseline until it was corrected [26] . Therefore, while the daily 
baseline repeatability provides a measure of consistency for 
orbits determined from multiday arcs, agreement with VLBI 
is a measure of overall system accuracy , which depends on a 
number of factors in addition to orbit accuracy. 

When different GPS antennas were used at the ends of a 
baseline, although daily repeatability was excellent, agreement 
with VLBI was worse by 1 to 7 cm, with no apparent depen- 
dence on baseline length. This was noticed in both the Novem- 
ber 1985 and June 1986 experiments. However, baselines with 
the same type of GPS antennas showed good agreement with 
VLBI. Since ephemeris errors tend to scale with baseline 
length, it was hypothesized that these discrepancies were due 
to local phenomena rather than orbital effects. Attention has 
been directed at phase center variations in the antennas, since 
the TI antennas are designed so that in operation the phase 
center variations nearly cancel out between sites that are not 
more than a few thousand km apart [6] . However, the Series-X 
antennas do not have the same phase center characteristics as 
the TI antennas, and the signatures resulting from the several- 
cm phase center variations that have been measured [27] 
could corrupt baseline measurements between unlike antennas. 
Therefore we qualify our high-precision results with the warn- 
ing that measurements between different GPS antennas may 
be much less reliable and may be affected by unpredictable 
effects. 

D. Baseline Repeatability Outside the 
Fiducial Network 

The baselines in North America are fairly well determined 
because they are either inside or near the fiducial network and 
because the current limited GPS constellation by design is op- 
timized for North America, especially the southwest United 
States. The formal errors from the filter are consistent with 
the results in Figs. 7, 8, and 9, predicting precision of 1 to 4 cm 
over these 2000-km baselines. To further test the robustness 
of the multiday-arc GPS orbits, we have determined baselines 
between Richmond and several sites in the Caribbean Sea 
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region occupied during the June 1986 experiment. Figure 10 
shows daily baseline repeatability for Richmond-Grand Turk 
(1049 km) and Richmond-Isabela (1582 km) determined 
from an 8-day orbit fit for June 2-10, 1986. Because the 
Caribbean sites are far to the southeast of both the fiducial 
network and the optimal region for the GPS constellation, the 
formal errors for these baselines are typically 2 to 7 cm, some- 
what worse than those of the North American baselines. De- 
spite the degraded geometry and reduced common visibility, 
the baselines to the Caribbean sites show precision of 1 to 4 
parts in 10 8 . Some of the Caribbean sites were equipped with 
WVRs; for these as well as the sites without WVRs, residual 
tropospheric corrections determined from the GPS data were 
critical in achieving these levels of baseline precision. Since the 
Caribbean sites are considerably more humid than most of the 
North American sites, various strategies for reducing system- 
atic errors due to uncertainties in the wet troposphere cor- 
rection are being studied with this data set. 

V. Conclusions 

It has been demonstrated that submeter GPS orbits can be 
determined using multiday arc solutions with the current GPS 
constellation subset visible for about 8 hours each day from 
North America. Submeter orbit accuracy was shown through 
orbit repeatability and orbit prediction. North American base- 
lines of 1000 to 2000 km in length can be estimated simul- 
taneously with the GPS orbits to an accuracy of better than 
1.5 parts in 10 8 (3 cm over a 2000-cm distance) with a daily 
precision of 2 parts in 10 8 or better. The most reliable baseline 
solutions are obtained using the same type of receivers and 
antennas at each end of the baseline. Baselines longer than 


1000 km between Florida and sites in the Caribbean region 
have also been determined with daily precision of 1 to 4 parts 
in 10 8 . The Caribbean sites are located well outside the fidu- 
cial tracking network and the region of optimal GPS common 
visibility, so these results further demonstrate the robustness 
of the multiday-arc GPS orbit solutions. 

Process noise models have been used in the orbit determina- 
tion filter to minimize systematic errors which can seriously 
affect ephemeris and baseline accuracy. These systematic 
effects include tropospheric delay fluctuations and small, un- 
modeled spacecraft accelerations. The process noise tropo- 
sphere models improved all orbit and baseline solutions, re- 
gardless of length of data arc. Tightly constrained process 
noise representation for part of the solar pressure model sig- 
nificantly improved baseline repeatability for arcs longer than 
1 week; however, an equally effective technique had fictitious 
thrusts estimated stochastically for each GPS satellite. Because 
of the limited ground visibility with the current constellation, 
it is not yet possible to determine whether the accelerations 
are genuinely related to solar radiation pressure or are due to 
other random or systematic forces acting on the spacecraft. 

This demonstration of several-cm accuracy over distances 
of a few thousand km, despite a limited ground tracking net- 
work and a constellation of only seven satellites, proves that 
GPS provides a very powerful relative positioning capability. 
It shows that GPS techniques have the intrinsic data strength 
and robustness needed for DSN high-precision applications, 
as well as low Earth orbiter tracking and crustal motion 
studies. 
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INDEPENDENT GPS ORBITS 

Fig. 3. Orbit repeatability for November 1985 
uses data from Nov. 13, 15, 18 for one solution 
and Nov. 14, 16, 19 for the other. Rms difference 
between the two solutions is computed over 
Nov. 17, during which no data was taken. 


Fig. 1. Locations of ground tracking sites used in the analysis of 
GPS data from the November 1985 and June 1986 experiments. 



Fig. 2. Daily baseline repeatability with 2-week orbit arcs 
(Nov. 13-24, 1985) for the Mojave-Fort Davis 1314-km baseline 
showing dramatic improvement when stochastic force parameters 
are estimated for GPS satellites. 
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Fig. 4. Rms orbit repeatability over a 6-hr interval on Nov. 17. Rms 
between independent arcs shown in Fig. 3 for altitude, cross-track, 
and down-track components have been averaged for all seven GPS 
satellites. 
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Fig. 5. Submeter orbit repeatability for GPS 6 and 8. Rms difference 
between two independent solutions is computed over a 24-hr 
interval on Nov. 17 when no measurements were used. 
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AND ARC 2 ORBITS 

Fig. 6. Orbit prediction test for Nov. 1 985. Orbits determined from Nov. 1 3-1 9 are mapped ahead, 
and the rms difference is computed with an independent solution determined from Nov. 20-24. 


2 


E 

llT 

o 

m 1 

DC 

LU 

LL 


Q 

Z 

2 o 

H 

co 

O 

a. 


-1 


Fig. 7. Difference between arc 1 predicted orbit and arc 2 orbit as 
shown in Fig. 6 for GPS 8. Rms is taken over a 6-hr interval. 
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Fig. 8. Orbit prediction difference similar to Fig. 7 except that here | 

the difference and rms are shown for a prediction interval of more * 

than 3 days. 






Fig. 10. Daily GPS baseline repeatability and agreement with VLBI 
for the Richmond-Haystack 2046-km baseline determined with a 
multiday orbit fit from Nov. 13-24, 1985. 


Fig. 9. Daily GPS baseline repeatability and agreement with VLBI 
for the Hat Creek-Fort Davis 1933-km baseline determined with a 
multiday orbit fit from Nov. 13-24, 1985. 
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Fig. 11. Daily GPS baseline repeatability for two baselines (1049 
and 1582 km) to the Caribbean region determined with multiday arc 
orbits from the June 1986 experiment. 



